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We calculate primordial abundances of light elements produced during big bang nucleosynthesis 
when the fine structure constant and/or the cosmic expansion rate take non-standard values. We 
compare them with the recent values of observed D, ''He and ^Li abundances, which show slight in- 
consistency among themselves in the standard big bang nucleosynthesis scenario. This inconsistency 
is not solved by considering either a varying fine structure constant or a non-standard expansion 
rate separately but solutions are found by their simultaneous existence. 
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Big Bang nucleosynthesis (BBN) theory calculates how 
much light elements are produced in the early universe. 
The standard BBN takes initial amount of baryons as 
only one input which is parametrized by baryon number 
density divided by photon number density: rj = nt/n-y. 
With ry ^ O(10~^°), the theory successfully predicts ob- 
served D, ^He and ^Li abundances extending over ten 
digits. 

However, the recent measurements of the primordial 
light elements abundances indicate that the success does 
not seem to be perfect 0, S S 13 ■ Such discrepancy 
between the observation and the theory is most likely as- 
cribed to the existence of unknown systematic errors in 
the observation. It is usually considered that systematic 
errors in D observation are smaller than those in "^He and 
^Li because D is observed in primordial objects, quasar 
absorption systems, but regression with respect to metal- 
licity is necessary to deduce primordial ^He and ''Li abun- 
dance. In addition, the observed D abundance is consis- 
tent with the baryon density from the cosmic microwave 
background (CMB) data which supports its robust- 
ness. Therefore, from this viewpoint, unexplored system- 
atic errors in both "^He and ^Li measurements solve the 
discrepancy [sif. 

From another viewpoint, the investigations correctly 
estimate the systematic errors in ^He and ^Li measure- 
ments so that the discrepancy is solved by non-standard 
physics. The recent studies on non-standard BBN in- 
clude non-standard expansion rate (number of neutrino 
species other than 3) |3, lepton asyinmetry and a 
varying fine structure constant, a [l^l- They can solve 
discrepancy between either D and ^He or D and ^Li but 
solution for three elements together is not obtained by 
their individual application j32| . 

In this paper, we show the current measurement of the 
three light elements is consistent without invoking further 
observational systematic errors if a is higher than to- 
day's value during BBN and the expansion rate is slower 
than the standard value. One might think that three 
parameters {rj, a and the expansion rate) necessarily ex- 
plain any three observations, but since the combination 
of non-standard expansion rate and lepton asymmetry 
(which has been investigated in Ref. |^) can only solve 
D-*He discrepancy, it is worth searching some combina- 
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FIG. 1: Standard BBN calculations of ^He, D and '^Li abun- 
dances as functions of rj are indicated by three curves whose 
width shows theoretical la uncertainty. The observational 
la uncertainties are expressed by the vertical extension of 
the boxes. They are drawn to overlap the theory curves so 
that their horizontal extension shows allowed range of rj. The 
larger box for *He is from Ref. [13 ] and the smaller from 
Ref. El 



tion to reconcile the three. Moreover, varying a and non- 
standard expansion of the universe may both appear from 
common models based on string theory which accommo- 
dates both a dynamical origin of coupling constants and 
unusual characteristics of spacetime such as extra dimen- 
sions. 

First of all, we summarize in Fig. ^ the current status 
of predicted and measured abundances of D, ^He and ^Li. 
Theoretical uncertainties are computed through Monte 
Carlo simulations using the values of Ref. based on 
the reaction rates of Ref. J^] . Measured values are taken 
from Refs. Q (Eq. d)) and [11 (Eq. ©) for "He, from 
Ref. [3 for D, and from Ref. [13 for ^Li: 

YiHcTO = 0.238 ±0.002 ±0.005, (1) 
Y4HciT = 0.2421 ±0.0021, (2) 
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FIG. 2; a-dependence of light element abundances. The cases 
for Aq/qo = (a — ao)/oio — 0, +0.05, —0.05 are drawn with 
solid, dashed, and dotted lines. 
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In Eq. the first uncertainty is statistical and the 
second one is systematic. Their root-mean-square, 
[(stat.)^ -I- (syst.)^]"'^/^, is the combined la error. For 
asymmetric errors, we adopt conservatively the larger one 
as Icr error (for ^Li, we divide the error in Eq. by 2 
to make it Icr). From the figure, we see ''He and ^Li are 
compatible with 77 w (2 ~ 4) x 10"^'^ but higher baryon 
density 77 « 6 x 10^^" is necessary for D. On performing 
analysis, due to the severer D-^Li discrepancy, we do 
not have a r] range to explain three elements abundances 
together with standard BBN at 99% confidence level (for 
either *He observations). We stress once again that such 
a discrepancy first requires a reassessment of systematic 
effects in the measurements of primordial abundances 
(especially that of ^Li), but below, assuming further sys- 
tematic errors are not found, wc investigate whether this 
discrepancy is solved by considering varying a and/or 
non-standard expansion rate. 

We calculate BBN abundances when a is different from 
the measured value at present, an ~ 1/137, as is de- 
scribed in our previous paper O [23. We have added 
the Ref. 's improvements concerning a-dependence of 
the nuclear reaction rates and binding energies. Fig. [21 
shows how abundances change when a is varied, repro- 
ducing the results of Ref. [13 . The dependence is un- 
derstood as follows. For ''He, since increasing a de- 



creases the neutron-proton mass difference (a proton is 
electrically charged so it becomes heavier than a neutron, 
which is electrically neutral), Am, the freeze-out ratio of 
neutron to proton increases and so does ^He abundance 
j2^j21j. Meanwhile, other light elements abundances are 
affected mainly by the change in the Coulomb barrier 
penetrability for the charged-particle induced nuclear re- 
action rates J,(3, ,22], which is the exponential factor in 
the following expression of the cross section (j{E) at en- 
ergy E, 



a{E) = cxp ( -27raZ,:Z, 

E 



(5) 



where S{E) is the astronomical S*- factor, fi is the reduced 
mass, and Zij are the atomic number of the colliding nu- 
clei. Since larger a suppresses the charged-particle reac- 
tion rates, the nucleosynthesis proceeds slower and this 
saves more D to be burned out. The same is true for 
T and more of it survives with higher a. This explains 
the ^Li increase for lower rj since ^Li is mainly produced 
by ^He(T,7)^Li. For higher 77, ^Li comes from the elec- 
tron capture of ^Be. ^Be is in turn produced through 
'*He('^He,7)^Be, which is strongly suppressed by higher 
a because of the large Coulomb barrier (Zi = Zj ~ 2 in 
Eq. © ). 

Especially, the dependence of "^He on a is derived from 
a number ratio of neutron to proton when their inter- 
change freezes out, that is, when the weak interaction be- 
comes comparable to the expansion rate of the universe. 
Since almost every neutron is synthesized into ^He, its 
mass fraction is approximately expressed by neutron and 
proton number density, n„ and Up at freeze-out (denoted 
by subscript "/") a-S 
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where freeze-out temperature Tf is about 0.7 MeV. Since 
Am is measured to be 1.293 MeV and its electromagnetic 
part is -0.76 MeV [H, 
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2.05 MeV. 
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Then AYif^^/Yifj^ « Aa/ao follows. 

To figure out whether varying a can solve the discrep- 
ancy between D and ^He and/or ^Li, we calculate as a 
function of 77 and a and search parameter space allowed 
by the observation of the light elements. The results are 
summarized in Fig. [3 We note that we take into ac- 
count the uncertainty in the present value of the electro- 
magnetic part of Am (which we neglected in Ref. [IB]). 
Since Ref. (231 have reported it to be less than 0.3 MeV, 
we regard this value to be 3a of a gaussian error pro- 
file and incorporate in our Monte Carlo simulation along 
with uncertainties in the reaction rates. This uncertainty 
does not exist when a = ao, at which theoretical errors 
become discontinuous and we have to split the calcu- 
lation like Fig. [3 
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FIG. 3: Contours show 95% confidence regions allowed by 
combinations of two elements observations, D and *He, D 
and ^Li, and *He and ^Li. As for combinations including 
*He, solid lines use *He data of Ref. and dotted lines use 
those of Ref. A region consistent with three elements 

abundances together is not found. We see D and ^Li are only 
reconciled by adopting a ^ ao but that makes ''He abundance 
inconsistent with the observation. 



In Fig. 131 as is expected from Figsn]and|21or has been 
demonstrated in Ref. , there is no region that explains 
measurements of the three elements together. Roughly 
speaking, ^He observation constrains in a-direction while 
D and ""Li constrain 77-direction because they are more 
sensitive to corresponding parameters (see Fig|2l. The 
^IIe-|-''Li contour lies around a — ao because they are al- 
ready consistent in the standard BBN. The D-|-^He con- 
tour hes in a < QfQ because, as can be seen in Fig. ^ 
^He is too much synthesized around the rj range deter- 
mined by D and cutting it down by decreasing a makes 
it consistent with the observation. At last, the D-f^Li 
contours exist in both a > ao and a < ao- The former 
corresponds to the region where D determines the 77-range 
(in which ^Li is over-synthesized) and the a-direction is 
constrained by required increase in a to bring down the 
^Li. Meanwhile, for the latter, since ^Li fix the 77-range, 
a < Qfo is necessary to decrease over-produced D. How- 
ever, such a ^ ao capable of reconciling D and ^Li either 




FIG. 4: A^^-dependence of light elements abundances. The 
cases for = 3, 4, 2 are drawn with solid, dashed, and 
dotted lines. 



over- or under-produce ^He. 

The story so far on BBN with varying a (and with the 
observations we adopt) is summarized as follows. It fails 
to explain the observed three light elements abundances 
because a required to make D and ''Li compatible creates 
too much or too small '^He. Since it is difficult to come up 
with any non-standard BBN other than varying a which 
reconciles D and ^Li, considering one which recovers "^He 
without violating the success of the varied a on D and 
^Li would be next best option. 

Our choice in this paper is non-standard expansion 
rate. It is usually treated and parameterized as effective 
number of neutrino species, N^, and we follow this con- 
vention. The standard BBN corresponds to iV^ = 3. The 
dependence of the light elements abundances on N^, is 
shown in Fig.^ which we explain briefly. "^He is again de- 
termined by Eq. ©. Since increasing N^, means increas- 
ing expansion rate, it raises freeze-out temperature and 
leaves more neutrons which synthesize into ^He. Faster 
expansion rate also makes nucleosynthesis less effective 
so more D is left unburned. So is T which fuses with "^He 
to be ^Li for lower rj. As for ''Be (which is the origin of 
^Li for higher 77), its destruction process ^Be(n,p)^Li is 
enhanced due to increased n. The reason for increase in 
n is same as D. 

We see from the figure that N^, does worse than a as for 
reconciling D and ^Li because barely alter them (es- 
pecially for higher 77) while changing the "^He abundance 
substantially. However, this turns out to be the advan- 
tage when combined with varying a because that is the 
very required property mentioned above. For a > ao, 
since ^He is overproduced in order to adjust D and ^Li, 
Ni, has to decrease. On the other hand, for a < ao, N^, 
has to increase. 

With this insight, we perform calculations similar 
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FIG. 5: Contours show 95% confidence regions allowed by 
combinations of all three elements observations, D, *He and 
^Li, for various A^,^. Solid lines use ^He data of Ref. 15] and 
dotted lines use those of Ref. . We do not find an allowed 
region for a < ao by varying A^^ . 

to Fig. O making N,j less than 3 for a > and more 
than 3 for a < a^. Naive expectation is that the contours 
including ^He seen in Fig. |3| would approach D+^Li con- 
tours and eventually they would merge to form allowed 
regions of three elements observations for both a > 
and a < ttQ. However, while solutions are obtained for 
a > ao as Fig. [3 there is no solution for a < a^. Notice 
that we look for solutions in the range 0.71 < a < 1.29, 
where the modification to the reaction rates caused by 
varying a is considered to be valid |1C| . The different be- 
havior stems from subtle effects of varying A''^ on D and 
'^Li. With Ni, < 3, from Figs. and ||| we see that the 
measured D and ''Li are more consistent (with smaller 
rj) than the standard BBN because D is predicted to 
be smaller (^Li is predicted to have larger abundance 
which means working in opposite way but N^- and 77- 
dependences are both larger for D, especially when the 
measured value of ^Li is around the trough of the theo- 
retical curve, so D is thought to be the decisive factor). 
On the contrary, since iV^ > 3 makes D larger, higher 77 
is more consistent. This requires very small a to decrease 



D, so small as to be outside the region with theoretical 
reliability. 

Fig.jSjshows our solutions to the inconsistency between 
the standard BBN and the measured primordial abun- 
dances of the light elements, Eqs. Q to The different 
^He measurements give similar results because the theo- 
retical uncertainty is comparable to the combined uncer- 
tainty in the observation of Eq. This mainly comes 
from the uncertainty in the electromagnetic part of the 
neutron-proton mass difference because the effect of un- 
certainties in the reaction rates on *He yield is negligible. 
The estimation for this error affects the size of the allowed 
regions but central values and qualitative features do not 
change. It is concluded that larger a and slower expan- 
sion rate (expressed by N^, < 3) solve the discrepancy 
between the standard theory and the observations. The 
solution with maximum N^, is found at about N,y « 1.16, 
77 « 4.7 X 10-^° and a/ao « 1.05. 

As we finish, we would like to make some comments. 
The first is on possible origins of lower-than-standard 
expansion rate (A'';^ < 3). We present here three possi- 
bilities: 1) non-thermal distribution of active neutrinos 
caused by low reheating temperature 0113; 2) negative 
dark radiation possibly exists in brane world scenarios 
; and 3) a varying (smaller at BBN) gravitational 
constant which is often found when the Ricci scalar is 
non-minimally coupled to a scalar field The 1) is 

limited to A^^ > but 2) and 3) can be any value as far as 
the total energy is positive. Their connection to varying 
a would be promising and quite interesting. The second 
is consistency with the CMB data. To our knowledge, 
there is no analysis of CMB data concerning simultane- 
ous change in a and N,y. Since full statistical treatment 
is beyond the scope of this paper, we just check the ef- 
fect of our solution on the first peak. An increase in a 
raises the first peak and a decrease in Ni, lowers it 
l^oj. It is reassuring that some cancellation is likely to 
take place but details remain to be worked out. Third, 
another solution should be found by considering lepton 
asymmetry instead of the non-standard expansion rate 
because its existence change *He abundance considerably 
while leaving D and ''Li almost unchanged. These issues 
are discussed in other places [30l |. 
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